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The global phase diagram of the spinless Falicov-Kimball model in d — 3 spatial dimensions is 
obtained by renormalization-group theory. This global phase diagram exhibits five distinct phases. 
Four of these phases are charge-ordered (CO) phases, in which the system forms two sublattices with 
different electron densities. The CO phases occur at and near half filling of the conduction electrons 
for the entire range of localized electron densities. The phase boundaries are second order, except 
for the intermediate and large interaction regimes, where a first-order phase boundary occurs in the 
central region of the phase diagram, resulting in phase coexistence at and near half filling of both 
localized and conduction electrons. These two-phase or three-phase coexistence regions are between 
different charge-ordered phases, between charge-ordered and disordered phases, and between dense 
and dilute disordered phases. The second-order phase boundaries terminate on the first-order phase 
transitions via critical endpoints and double critical endpoints. The first-order phase boundary is 
delimited by critical points. The cross-sections of the global phase diagram with respect to the 
chemical potentials and densities of the localized and conduction electrons, at all representative 
interactions strengths, hopping strengths, and temperatures, are calculated and exhibit ten distinct 
topologies. 

PACS numbers: 71.10.Hf, 05.30.Fk, 64.60.De, 71.10.Fd 



I. INTRODUCTION 

The Falicov-Kimball model (FKM) was first proposed 
by L. M. Falicov and Kimball [l| to analyze the ther- 
modynamics of semiconductor-metal transitions in SmBg 
and transition- metal oxides The model incor- 

porates two types of electrons: one type can undergo 
hopping between sites and the other type cannot hop, 
thereby being localized at the sites. Thus, in its introduc- 
tion, FKM described the Coulomb interaction between 
mobile d band electrons and localized / band electrons. 
There have been a multitude of subsequent physical in- 
terpretations based on this interaction, including that of 
localized ions attractively interacting with mobile elec- 
trons, which yields crystalline formation [y, |7|. Another 
physical interpretation of the model is as a binary alloy, 
in which the localized degree of freedom reflects A or B 
atom occupation [H HI . In this paper we employ the orig- 
inal language, with d and / electrons as conduction and 
localized electrons with a repulsive interaction between 
them. 

Since there is no interacting spin degree of freedom 
in the Hamiltonian, the model is traditionally studied 
in the spinless case, commonly referred as the spinless 
FKM (SFKM) and which is in fact a special case of the 
Hubbard model in which one type of spin {e.g., spin-up) 
cannot hop [Tof. In spite of its simplicity, this model is 
able to describe many physical phenomena in rare-earth 
and transition metal compounds, such as metal transi- 
tions, charge ordering, etc. 

Beyond the introduction of the spin degree of free- 



dom for both electrons [TTl - [27l |. there also exist many 
extensions of the original model. The most widely stud- 
ied extensions include multiband hybridizatio n I28l435| . 
/ — / hopping [3^53, correlated hopping [41H45I ]. non- 
bipartite lattices|4f 4471 . hard-core bosonic particles [13] , 
magnetic fields [IS, 12414271 l47l |48| . and next- nearest- 
neighbor hopping 491. Exhaustive reviews are avail- 
able in Refs. [50l453l |. The wider physical applica- 
tion of both the basic FKM and its extended versions 
have aimed at explaining valence transitions [IH, [l8l - [20| , 
metal-insulator transitions [lU, [2l| - [23l |54| , mixed valence 
phenomena [55|. Raman scattering [56| . colossal mag- 
netoresistance [2414271 , electronic ferroelectricity [3J, 1371 - 
Hi| , and phase separation [U, IH 5J HHHI • 

After the initial works on the FKM pHll, the litera- 
ture had to wait 14 years for the celebrated first rigorous 
results. Two independent studies, by Kennedy and Lieb 
[1, 0] and by Brandt and Schmidt [6(J Hl|, proved for 
dimensions d > 2 that, at low temperatures, FKM has 
long-range charge order with the formation of two sublat- 
tices. Various methods have been used in the study of the 
FKM. In most of these studies, either the d^oo infinite- 
dimensional limit or d = 1,2 low-dimensional cases have 
been investigated. Studies include limiting cases such 
as ground-state analysis or the large interaction limit. 
Renormalization-group theory [62| offers fully physical 
and fairly easy techniques to yield global phase diagrams 
and other physical phenomena. 

This non-trivial nature of SFKM motivated us to de- 
termine the global phase diagram of the model, which re- 
sulted in a richly complex phase diagram involving charge 
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FIG. 1: Evolution of the cross-sections of the global phase diagram under t increase for \U\ = 1, in terms of the chemical 
potentials (upper panels) and densities (lower panels) of the localized and conduction electrons. In phase subscripts throughout 
this paper, the first and second subscripts respectively describe localized and conduction electron densities, as dilute (d) or 
dense (D). The dotted and thick full lines are respectively first- and second-order phase transitions. Phase separation, i.e., 
phase coexistence occurs inside the dotted boundaries, as identified appropriately but not repeatedly in the figure. The dashed 
lines are not phase transitions, but smooth changes between the different density regions of the disordered (S) phase. The 
charge-ordered phases are denoted by CO. The charge-ordered phases occur as strips near the half filling of the conduction 
electrons. Phase separation occurs near the simultaneous half filling of both the localized and conduction electrons. The four 
rounded coexistence regions, two of which disappear as t is increased, are two-phase coexistence regions between the disordered 
S phases that are distinguished by electron densities. The narrow triangular regions are three-phase coexistence regions between 
the S phases. The second-order transition lines bounding the charge-ordered CO phases terminate at critical endpoints on the 
coexistence regions. These endpoints, as t is increased, move past each other, as detailed in Sec. IV below, leading to coexistence 
regions between the different charge-ordered phases and between the charge-ordered and disordered phases. 



ordering and phase coexistence, as exemplified in Fig. 1. 

We use the general method for arbitrary dimensional 
quantum systems developed by A. Falicov and Berker [631 ] 
to obtain the global phase diagram of the SFKM in d = 3, 
in terms of both the chemical potentials and the densities 
of the two types of electrons, for all temperatures. The 
outline of this paper is as follows: In Sec. II, we introduce 
the SFKM and, in Sec. Ill, we present the method [63|. 
Calculated phase diagrams are presented in Sec. IV, for 
the non-hopping (i = 0) classical submodel and for the 
hopping (i ^ 0) quantum regimes of small, intermediate, 
and large \U\. We conclude the paper in Sec.V. 

II. SPINLESS FALICOV-KIMBALL MODEL 

The SFKM is defined by the Hamiltonian 

-p% = t (c]cj + ct Ci ) + U ^2 n lWi 

Hi) < (1) 

+ no m + v Wj , 

i i 

where ft — 1/fegT and (ij) denotes that the sum runs 
over all nearest-neighbor pairs of sites. Note that, as 
in all renormalization-group studies, the Hamiltonian 



has absorbed the inverse temperature. The dimension- 
less hopping strength t can therefore be used as the in- 
verse temperature. Here cj and ci are respectively cre- 
ation and annihilation operators for the conduction elec- 
trons at lattice site i, obeying the anticommutation rules 
{c l; Cj} = {c\, c]} = and {cj, Cj} = while rn = c\a 
and Wi are electron number operators for conduction and 
localized electrons respectively. The operator Wi takes 
the values 1 or 0, for site i being respectively occupied 
or unoccupied by a localized electron. The particles are 
fermions, so that the Pauli exclusion principle forbids 
the occupation of a given site by more than one localized 
electron or by more than one conduction electron. 

The first term of the Hamiltonian is the kinetic energy 
term, responsible for the quantum nature of the model. 
The system being invariant under sign change of t (via a 
phase change of the local basis states in one sublattice), 
only positive t values are considered. The second term 
is the screened on-site Coulomb interaction between lo- 
calized and conduction electrons, with positive and neg- 
ative U$ values corresponding to attractive and repulsive 
interactions. We consider only the repulsive case, since 
the attractive case can be connected to the repulsive one 
by the particle-hole symmetry possessed by either type of 
electrons. Particle- hole symmetries are achieved by the 
transformations of Wi — >■ 1 — lOj for the localized electrons 
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and c. 



KjcJ for the conduction electrons, 



where, for a bipartite lattice, Ki = 1 for one sublatticc 
and Ki = — 1 for the other [8, 53]. The last two terms of 
the Hamiltonian are the chemical potential terms with v$ 
and /io being the chemical potential for a localized and 
conduction electron. 

In order to carry out a renormalization-group trans- 
formation easily, we trivially rearrange the Hamiltonian 
given in Eq.© into the equivalent form of 

—f3Ti = ^2 t (clcj + cjcjj + U (mwi + rijWj) 



+ [l (m + Tij) + V (Wi + Wj) 



(2) 



m 

where, for a d- dimensional hypercubic lattice, U = 
Uo/2d, [i = fMy/2d, v = i/o/2d, and is the two-site 

Hamiltonian involving only nearest-neighbor sites i and 

3- 



III. RENORMALIZATION-GROUP THEORY 



Suzuki- Takano Method in d = 1 



In d = 1 the Hamiltonian in Eq.© is 
-/3H = J2 [-PHi,»i] ■ 



(3) 



The renormalization-group procedure traces out half of 
the degrees of freedom in the partition function [rJil [65j j , 



Tr odd e-^ =Tr odd eSJ-^,* +1 ] 



-mi,i+i] 



(4) 



n« 



;£° dd [-/i'«--l, I + l] = e -/3'«'. 



Here and throughout this paper primes are used for the 
renormalized system. Thus, as an approximation, the 
non-commutativity of the operators beyond three con- 
secutive sites is ignored at each successive length scale, 
in the two steps indicated by ~ in the above equation. 
Earlier studies [631473 ] have established the quantitative 
validity of this procedure. 

The above transformation is algebraically summarized 

in 



where i,j,k are three successive sites. The opera- 
tor —p'H'i k acts on two-site states, while the operator 
~(3TLij — PT~Cj,k acts on three-site states. Thus we can 
rewrite Eq.© in matrix form as 



[UiV k \e 



Sj v k \e-^~^ \Ui sj v k ), (6) 



where state variables U£, vt, si, ut, and ve can be one of 
the four possible single-site \wi,nz) states at each site £, 
namely one of |00), |01), |10>, and |11). Eq.© indicates 
that the unrenormalized 64 x 64 matrix on the right-hand 
side is contracted into the renormalized 16 x 16 matrix on 
the left-hand side. We use two-site basis states, {|</> P }}, 
and three-site basis states, { | ) } ^ m order to block- 
diagonalize the matrices in Eq. © . These basis states are 
the eigenstates of total localized and conduction electron 
numbers. The set of {|0 P )} and {|^ g )} are given in Tables 
U and [H] respectively. The corresponding block-diagonal 
Hamiltonian matrices are given in Appendices A and B. 
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TABLE I: The two-site basis states that appear in Eq.lO, 
in the form \wiTii, WjUj). The total localized and conduction 
electron numbers w and n, the eigenvalue u of the operator 
Tij defined after Eq.© are indicated. |0g_i2) are respectively 
obtained from |05-s) by the action of Tij, while the corre- 
sponding Hamiltonian matrix elements are multiplied by the 
u values of the states. 



(5) 



With these basis states, Eq.© can be rewritten as 

(0p|e _/3 Hi > k \<j)p) = ^ ^2(^p\UiV k )(UiSjVk\tpq)- 
q,q 

u.v,s 

Once written in the basis states {|0 P )}, the block- 
diagonal renormalized matrix has 13 independent ele- 
ments, which means that renormalization-group trans- 
formation of the Hamiltonian generates 9 more interac- 
tion constants apart from t, U , fi, and v. In this 13- 
dimensional interaction space, the form of the Hamilto- 
nian stays closed under renormalization-group transfer- 
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mations. This Hamiltonian is 



TABLE II: The three-site basis states that appear in Eq.©, 
in the form \wini,Wjrij,Wknk). The total localized and con- 
duction electron numbers w and n, the eigenvalue u of the 
operator Tik defined after Eq.(j8j are indicated. {|f/>ii+3a;)}, 
x = 0,1, ... ,15, are respectively obtained from {1^)9+32,)} by 
the action of Tik , while the corresponding Hamiltonian matrix 
elements are multiplied by the u values of the states. 



-PT-Li,j = t (c\cj + ctcjj + U (riiWi + rijWj) 

+ /i (rij + rij) + v (wi + Wj) + Jriirij 
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+ V n n i n j (wi +Wj) + V w (n t + nj)w t Wj 



(8) 



QTijWiWj + RTij(wi + wj 



G, 



where T t j is a local operator that switches the conduc- 
tion electron states of sites i and j: Tij\wini,Wjnj) — 
u\winj,Wj7ii) with u = 1 for rii + rij < 2 and u 

I otherwise. When To- is applied, further below, to 
three consecutive sites i, j, k, Tik\wini,Wjrij,Wknk) = 
u\wink,Wjrij , Wkfii) with u = 1 for ni + rij + rik < 2 and 
u = — 1 otherwise. 



To extract the renormalization-group recursion re- 
lations, we consider the matrix elements 7 p .p = 

(0p|e" /3 Wi < fc |0p)- With7 9i9 = 75 iS , 710,10 = 76,6, 711,11 = 
77,7, and 712,12 = 78,8, 12 out of 16 diagonal elements are 
independent and, with 7 10) n = 7n )10 = -77^ = -76,7, 
only one of the 4 off-diagonal elements is independent, 
summing up to 13 independent matrix elements. Thus 
we obtain the renormalizcd interaction constants in terms 
of {7}, defining 7 P = 7 PiP for the diagonal elements and 
70 = 76,7 for the only independent off-diagonal element: 
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The matrix elements {7} of the exponentiated renormal- 
ized Hamiltonian are connected, by Ecj.©, to the matrix 
elements, rj q ^ = (ijj q \e~ l3 ' Hi ' i ~P' H -i- k \ipq) of the exponenti- 
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TABLE III: Interaction constants K a , runaway coefficients K' a /K a , and expectation values M a — (K a ), at the phase sinks. 
Here, K a are used as abbreviations for the conjugate operators for interaction constants K a , e.g., (t) = {c\cj + cjcj), (U) = 
{riiWi + iijWj), etc. The non-zero hopping expectation value is —a — —0.629050. In the subscripts in the first columns, the left 
and right entries refer to the localized and conduction electrons, respectively, as dilute (d) or dense (D). 



Phase 


Boundary Interaction constants K a at the boundary fixed points 


Relevant eigenvalue 


boundary 


type 


t 


U n v J K L P V„ V w Q 


R 




exponent yi 


CO dD /CO D d 


1st 
order 


oo 


— oo oo — oo oo oo oo — oo — oo oo oo 


— oo 


2/z - 2u + J - K 
-Q-2R = Q 


3 


COdd/^dd 


2nd 
order 


oo 


oo — oo — oo oo — oo oo — oo — oo oo oo 


— oo 


t + n = 1.744253 


0.273873 


COdD/^dD 


2nd 
order 


oo 


oo oo — oo oo — oo oo — oo oo oo oo 


— oo 


t-n-J= 1.744253 


0.273873 


COod/fod 


2nd 
order 


oo 


oo — oo oo oo — oo oo — oo — oo oo oo 


— oo 


t + U + n + P + V w 
= 1.744253 


0.273873 


COdd/5dd 


2nd 
order 


oo 


oo oo oo oo — oo oo oo — oo oo oo 


oo 


t-U-fi-J-L-P 
-2V„ -V W + 2Q + 4R = 


0.273873 


COdd/CO d D 


2nd 
order 


oo 


oo — oo — oo oo — oo oo — oo — oo oo oo 


— oo 


2(j, + J = 


1.420396 



CODd/CODD 2nd oo oo — oo oo oo — oo oo — oo oo oo oo — oo 2U + 2fj, + J + L + 2P 1.420396 
order +2V n + 2V W - 2Q - 4R = 



TABLE IV: Interaction constants K a and relevant eigenvalue exponents j/i at the phase boundary fixed points. For first-order 
phase transitions, j/i = d = 3. 

ated unrenormalized Hamiltonian, The matrix elements r\ q ^ can be obtained in terms of 

70 =»?12,18 + f/21,27 + 7/36,42 + ??45,51, 

71 =»7l +V3 + Via + Vl6, 

72 =m + V7 + »?13 + V28, 

73 =m + V5 + VW + V22, 
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FIG. 2: Renormalization-group flow basins of the t = clas- 
sical submodel, in the chemical potentials (upper panel) and 
densities (lower panel) of the localized and conduction elec- 
trons. In phase subscripts throughout this paper, the first and 
second subscripts respectively describe localized and conduc- 
tion electron densities, as dilute (d) or dense (D). The dashed 
lines are not phase transitions, but smooth changes between 
the four different density regions of the disordered (8) phase. 



the unrenormalizcd interactions via exponentiating the 
unrenormalized Hamiltonian matrix whose elements are 
given in Appendix B. 



B. Renormalization-Group Transformation in d > 1 



Equations ([9]) and (fT0|) , together with Appendix B, 
constitute the renormalization-group recursion relations 
for d = 1, in the form K' = R(K), where K = 



(t,U,fi,v,J,K,L,P,V n ,V w ,Q,R,G). To generalize to 
higher dimension d > 1, we use the Migdal-Kadanoff 
procedure [zEIzl], 



K' 



,d-i 



R(K), 



11) 



where b = 2 is the rescaling factor and R is the 
renormalization-group transformation in d = 1 for the 
interaction constants vector K . This pro cedure is exact 
for d-dimensional hierarchical lattices [77H79| and a very 
good approximation for hypercubic lattices for obtaining 
complex phase diagrams. 

Each phase in the phase diagram has its own (stable) 
fixed point (s), which is called a phase sink (Table ITTTT) . 
All points within a phase flow to the sink(s) of that phase 
under successive renormalization-group transformations. 
Phase boundaries also have their own (unstable) fixed 
points (Table IIV[) . where the relevant exponent analy- 
sis gives the order of the phase transition. Thus, the 
repartition of the renormalization-group flows determine 
the phase diagram in thermodynamic-field space. Matrix 
multiplications, along the renormalization-group trajec- 
tory, with the derivative matrix of the recursion relations 
relate the expectation values at the starting point of the 
trajectory to the expectations values at the phase sink. 
The latter are determined (Table III) by the left eigen- 
vector, with eigenvalue b d , of the recursion matrix at the 
sink, where b = 2 is the length-rescaling factor of the 
renormalization-group transformation. When the expec- 
tation values are thus calculated for the points of the 
phase boundary, the phase diagram in density space is 
determined. HI El 



IV. GLOBAL PHASE DIAGRAM OF SFKM 

The global phase diagram of SFKM is calculated, as 
described above, for the whole range of the interac- 
tions (i, U, v, /i). The global phase diagram is thus 
4-dimcnsional. 1 jt can be taken as the temperature vari- 
able. We present the calculated global phase diagram 
in four subsections: The first subsection gives the t = 
classical submodel. The other subsections are devoted to 
small, intermediate, and large values of the interaction 
\U\. We present constant t/\U\ cross sections in terms 
of the localized and conduction electron chemical poten- 
tials v/\U\ and /i/|J7|and in terms of the localized and 
conduction electron densities (u)i) and (rii). 

A. The Classical Submodel t = 

Setting the quantum effect to zero, t = 0, yields 
the classical submodel, closed under the renormalization- 
group flows. The global flow basins in v/\U\ and (i/\U\ 
are the same for all U, given in Fig. [2] There exist four 
regions of a disordered phase within this submodel, which 
are localized-dilute-conduction-dilute, localized-dilute- 
conduction-dense, localized-dense-conduction-dilute, and 
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FIG. 3: Constant t/\U\ cross-sections of the phase diagram for interaction \U\ = 0.1, in terms of the chemical potentials (upper 
panels) and densities (lower panels) of the localized and conduction electrons. In phase subscripts throughout this paper, the 
first and second subscripts respectively describe localized and conduction electron densities, as dilute (d) or dense (D). The full 
lines are second-order phase transitions. The dashed lines are not phase transitions, but smooth changes between the different 
density regions of the disordered (8) phase. The charge-ordered phases are denoted by CO. Details are shown in Fig. 4. Thus, 
for low values of the interaction, all phase boundaries are second order and there is no phase coexistence. 



localized-dense-conduction-dense regions, denoted by 
cW, <JdD, fed, and <jdd- [In phase subscripts through- 
out this paper, the first and second subscripts respec- 
tively describe localized and conduction electron densi- 
ties, as dilute (d) or dense (D).] In the renormalization- 
group flows, each S region is the basin of attraction of 
its own sink. The dashed lines between the different re- 
gions are not phase boundaries, but smooth transitions 
(such as the supercritical liquid - gas or up-magnetized 
- down-magnetized transitions), which are controlled by 
zero-coupling null fixed points. [82j 

It should be noted that the Suzuki- Takano and Migdal- 



Kadanoff methods are actually exact for this classical 
submodel, and yield exactly the same picture as obtained 




B. The Small \U\ Regime 

In this subsection, we present our results for \U\ = 0.1, 
representative of the weak-interaction regime. The t = 
phase diagram of Fig. 2 evolves under the introduction 
of quantum effects via a non-zero hopping strength t. It 
should be noted that increasing the dimensionless Hamil- 
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FIG. 4: Zoomed portion of Fig. 3, for the \U\ = 0.1, t/|[/| = 
10 phase diagram. 



tonian parameter t is equivalent to reducing temperature, 
as in all renormalization-group studies. The first effect 
is the decrease and elimination (left panels of Fig. 3) 
of the (smooth) passage between the fed and Sdu re- 
gions. With this elimination, all four regions meet at 
v/\U\ = n/\U\ = 0.5 and {w t ) = (n,) = 0.5, the half 
filling of both localized and conduction electrons. With 
increasing t (equivalent to decreasing temperature) , four 
new, charge-ordered (CO) phases emerge at t ~ 0.6. The 
CO phases occur at and near half filling of conduction 
electrons for the entire range of localized electron densi- 
ties. The CO phases grow with increasing t (decreasing 
temperature) until saturation at high t (right panels of 
Fig. 3). 

All of the new CO phases have non-zero hopping den- 
sity (c\cj +CjCi) = — a = —0.629050 at their phase sinks. 
The expectation values at the sinks are evaluated as the 
left eigenvector of the recursion matrix with eigenvalue b d 
[80| . In the CO phases, the hopping strength t diverges 
to infinity under repeated renormalization-group trans- 
formations (whereas in the S phases, t vanishes under 
repeated renormalization-group transformations). The 
localized electron density is (wi + Wj) = at the sinks 
of COdd and COdD, while (w; + Wj) = 2 at the sinks of 
CODd and COdd, which throughout the corresponding 
phases calculationally translates [8(3] as low (d) and high 
(D) localized electron densities, respectively. Recall that 
on phase labels (CO and 6) throughout this paper, the 
first and second subscripts respectively describe localized 
and conduction electron densities. 

The conduction electron density is (rij + nj) — a — 
0.629050 at the sinks of COdd and COod, while (rij + 
rij) = 2 — a — 1.370950 at the sinks of COdD and COdd- 
The nearest-neighbor conduction electron number cor- 
relation is {uifij) — at the sinks of COdd and COod, 
while {riiUj) = 1 — a = 0.370950 at the sinks of COdD and 
COdd- Consequently, for conduction electrons, if a given 
site is occupied, its nearest-neighbor site is empty at the 
sinks of COdd and COod- The COdD and COdd phases 
are connected to the COdd and COod phases by particle- 
hole interchange on the conduction electrons. Thus, in 



the CO phases, the lattice can be divided into two sub- 
lattices with different electron densities. The behavior at 
the CO sinks therefore indicates charge ordered phases 
at finite temperatures, as also previously seen in ground- 
state studies [6l|, [84j, [85J . Note that this charge ordering 
is a purely quantum mechanical effect caused by hopping, 
since the SFKM Hamiltonian [Eq.([T])] studied here does 
not contain an interaction between electrons at different 
sites. 

In the small \U\ regime, all phase boundaries around 
the CO phases are second order. As seen in the expanded 
Fig. 4, all four CO phases and all four regions of the 5 
phase (as narrow slivers) meet at v/\U\ = fJ-/\U\ = 0.5 
and (wi) = (rii) = 0.5, half- filling point of both localized 
and conduction electrons. All characteristics of the sinks 
and boundary fixed points are given in Tables Unl and HVl 



C. The Intermediate \U\ Regime 



In this subsection, the phase diagram for \U\ = 1, rep- 
resentative of the intermediate- interaction regime, is pre- 
sented. Fig. Ogives constant t/\U\ cross-sections. First- 
order phase boundaries appear in the central region of 
the phase diagram, at and near the half filling of both 
localized and conduction electrons. 

For low values of t (left panels of Fig. 5), equivalent 
to high temperatures, two first-order phase boundaries, 
bounded by four critical points C, pinch at a quadruple 
point Q. In the (left-lower) density-density phase dia- 
gram, four phase separation (coexistence) regions mark 
the first-order phase transitions. Inside these regions, co- 
existence (phase separation) occurs between the phases 
on each side of these regions, as indicated on the figure. 
The tie line of the quadruple point is shown as a thin 
straight line. All four S phases coexist (phase separate) 
on this line. 

As t increases (temperature decreases), the four 
charge-ordered CO phases appear again at t ~ 0.6, as 
seen in the leftmost panels of Fig. 1. The CO phases 
again occur at and near half filling of conduction elec- 
trons for the entire range of localized electron densities. 
In the right panels of Fig. 5, the second-order transition 
lines bounding the CO phases terminate at two critical 
endpoints E [821 ] and two double critical endpoints E2 on 
the first-order line in the central region (zoomed in Fig. 
6). Thus, first-order transitions and phase separation oc- 
cur between the pairs of fed and <5dd,fed and COdd,COod 
and COdD, COdd and (5dD, fe>D and <5dD phases, as in- 
dicated on Fig. 6, at and near the half filling of both 
localized and conduction electrons. 

The evolution of the phase diagrams between right and 
left panels of Fig. 5 are shown in Fig. 1. 
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FIG. 5: Constant t/\U\ cross-sections of the phase diagram for interaction \U\ = 1, in terms of the chemical potentials (upper 
panels) and densities (lower panels) of the localized and conduction electrons. The dotted and thick full lines are respectively 
first- and second-order phase transitions. Phase separation, i.e., phase coexistence occurs inside the dotted boundaries, as 
identified in the figure. The details of the coexistence region in the lower-right panel are given in Fig. 6. The quadruple point 
Q tie line is shown as the thin straight line. The dashed lines are not phase transitions, but smooth changes between the 
different density regions of the disordered (5) phase. 



D. The Large \U\ Regime 

The evolution of the global phase diagram, as the in- 
teraction strength is increased, is seen in the phase di- 
agrams in Fig. 7. The CO phases emerge again at 
t ~ 0.6. With increasing t (decreasing temperature), the 
CO phases grow, until saturation seen in Fig. 7. The 
topology of the phase diagram with five phases stays the 
same for all t > 0.6. 

The constant t/\U\ cross-sections of the phase diagram 
are given in Fig. 7. For U = 1.5, the double crit- 
ical endpoints E2 have split into pairs of simple criti- 



cal endpoints E, resulting in six separate critical end- 
points. For U = 1.845628, the inner two critical end- 
points have merged into a single double critical end- 
point. For U = 10, the double critical endpoint has 
split into two critical endpoints and the critical lines in 
the low-density and high-density localized electrons re- 
gions have disconnected from each other. In this strong 
interaction limit, the homogenous (non-phase-separated) 
charge-ordered phases occur again at and near half filling 
of conduction electrons, but at the low- or high-density 
limit of the localized electrons. Away from these limits, 
the charge-ordered phases occur in coexistence (phase- 
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separated from) the disordered phases. At and near the 
half filling of both localized and conduction electrons, 
the coexistence of the disordered phases fed and 5^ oc- 
curs. Two sets of three critical lines terminate in separate 
endpoints, as seen in the zoomed Fig. 9. In this case, a 
characteristic shape of the density phase diagrams, which 
we dub chimaera coexistence, emerges. In the chimaera 
phase diagram, coexistence can be found for essentially 
the entire range of conduction electrons densities or for 
most of the range of localized electron densities. In the 
upper-right chemical-potential panel of Fig. 7, the first- 
order phase boundary exhibits, from the left to right, a 
sequence of the maximum, minimum, maximum, min- 
imum points; the four corresponding tie lines are also 
shown in the lower-right density panel. These tie lines 
abut, on one end, very near maxima and minima of the 
lower and upper branches of the coexistence boundaries, 
thereby underpinning the distinctive chimaera topology. 



For the t = classical submodel, we have obtained 
disordered (S) regions, dilute and dense separately for 
localized and conduction electrons, but no phase tran- 
sition between them. The repartition of these regions, 
delimited by renormalization-group flows, quantitatively 
stays the same for the whole \U\ range and is exactly 
as obtained in Ref. [Hj]. For the whole \U\ range and 
< t < 0.6, the classical submodel phase diagram is per- 
turbed in such a way that regions S^d and feo intercede 
between regions S^b and fed, resulting in the shrinking 
and disappearing of the S^b to fed passage. 

All d regions have vanishing hopping density at their 
corresponding sinks. For the whole \U\ range, upon in- 
creasing t (lowering temperature), at t ~ 0.6 four new 
phases (CO) emerge with non-zero hopping density of 
—a = —0.629050 at their sinks. These CO phases are 
also either dilute or dense, separately, in the localized and 
conduction electrons (COdd, COdD, COod, and COdd) 
and are all charge ordered in the conduction electrons, a 
wholly quantum mechanical effect. In these CO phases 
the bipartite lattice is divided into two sublattices of al- 
ternating electron density. The CO phases occur at or 
near the half filling of conduction electrons. The phase 
diagrams with all five phases for t > 0.6 exhibit differ- 
ent topologies, for the small, intermediate, and large \U\ 
regimes: 

For the small \U\ (weak-interaction) regime, all phase 
boundaries are second order. All five phases meet at 
v/\U\ = n/\U\ = 0.5 and (w,) = (m) = 0.5, the half- 
filling point of both localized and conduction electrons. 

For the intermediate \U\ (intermediate- interaction) 
regime, a first-order phase boundary emerges in the cen- 
tral region of the phase diagram. This first-order bound- 
ary is centered at v/\U\ = n/\U\ = 0.5 and is bounded 
by two critical points C. The second-order lines bound- 
ing the CO phases terminate at critical endpoints E and 
double critical endpoints E2 on the first-order boundary. 
Due to this first-order phase transition at and near the 
half filling of both localized and conduction electrons, a 
rich variety of phase separation (phase coexistence) oc- 
curs, as indicated on Figs. 1,5,6,7. 

For the large \U\ (strong-interaction) regime, as \U\ is 
increased, the critical endpoints pass through each other 
by merging and unmerging as double critical endpoints. 
For large \U\, the COod and COdd phases are detached 
from the COdd and COdD phases, forming two separate 
bundles, at high- and low-densities of localized electrons 
respectively. First-order transitions occur between the 
variously dense and dilute S. The global phase diagram 
underpinning all of these cross-sections is decidedly quite 
complex. 



CONCLUSION 
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APPENDIX A: BLOCK-DIAGONAL 
RENORMALIZED HAMILTONIAN 

The matrix elements of the block-diagonal renormal- 
ized 2-site Hamiltonian in the {|</> p )} basis are given in 
Eq.lJE), where (</> p \ - P'W^p) = e p + G for the 12 in- 
dependent diagonal elements and (^l — P'H'i tl^) = e o 
for the only independent off-diagonal element: 

e% = 0, 62 =t' + /jf, 63 = —t' + n', e 4 = 2// + j', 
e 5 = v' + R', e e = t' + U'/2 + fi' + v' + P'/2 + R', 
e 7 = ~t' + U'/2 + // + v' + P'/2 - R', 
6 8 = U' + 2/i' + v' + J' + P' + - R', 
6 13 = 2v' + K' + Q' + 2R', 
6 14 =t' + U' + fi' + 2v' + K' + P' + Vl + Q' + 2R', 
ei5 = + U' + p' + 2i/ + K' + P' + VI - Q' - 2R', 
ei6 = 2(U'+/j l '+i/)+J'+K'+L'+2(P'+V n +Vl)-Q'-2R', 
6 = (W - P')/2. 

(12) 



The matrix elements for the states connected by the ex- 
change of the outer conduction electrons are obtained 
by multiplication with the eigenvalues u of T^. The 
matrix elements rj q ,q that enter the recursion relations 
via Eq. (fT0| arc obtained by exponentiating the block- 
diagonal Hamiltonian given here. 



APPENDIX B: BLOCK-DIAGONAL 
UNRENORMALIZED HAMILTONIAN 



The matrix elements of the block-diagonal unrenormal- 
ized 3-site Hamiltonian in the {iV'g)} basis are given in 
Eq.CESl), where (ip q \- j3Hij- PU jtk \ip q } = e q + 2G for the 
diagonal elements and (-0^ | — /3Hij — /3'Hj j k\ipq) — £ q ,q for 
the off-diagonal elements: 

£l = 0, 62 = £3 = £4 = £6/2 = /i, £5 = £7 = 2/1 + J, 

£s = 3/! + 2J, £9 = £34/2 = v + R, Ew = v + 2R, 
£12 = eis = U/2 + n + v + R/2, 
£13 = £19 = ^ + v + P + R, £15 = p + u + P, 

£l6 = £49/2 = U + I-J. + V, 

£21 = U/2 + 2fi + is/2 + J + P+ (V„ - R)/2, 
£22 = £28 = U + 2/J, + V + J + P + V n - R, 

£24 = U + 2/i + v, £ 25 = 2/i + v + 2P, 
£27 = U/2 + 2fi + i/ + J + P+ (V n - R)/2, 
£30 = U + 3fj, + v + 2 J + P + V„ - R, 
e 3 i = U + 3/i + v + 2(J + P + V n - R), 
e 3a = 2v + K + Q + 3i?, 
£ 36 = £42 = U/2 + n + 2u + K + P + (y w +Q + 3i?)/2, 

£37 = £43 = U + /i + 2v + R, 
e 39 = U + fi + 2v + K + P + V w , £ 40 = fi + 2(v + P), 
e 45 = eB1 = 3U/2 + 2(/x + u) + J + K + L/2 + 2P 

+ 3(F„ + V w )/2 - (Q + 3R)/2, 
£46 = £52 = U + 2(/x + v) + J + 2P + V n - R, 
£ 48 = U + 2(/z + v) + K + 2P + V w , 
£54 = 2U+3fi+2(is+J)+K+L+3(P+V n )+2V w -Q-3R, 
£55 = 2U + 3/i + 2[v + J + P + V n - R), 
£57 = iv + 2(K + Q) + 4R, 
e 5 8 = e 6 Q = U + (x + 3v + 2K + P + V w +Q + 2R, 
£ KQ = U 4- 11 + 3i/ + 2(K + P + V..,). 
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